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ABSTRACT 


The rotational spectral lines of C-C 3 H 2 and two kinds of the 13 C isotopic 
species, c- 13 CCCH 2 (C' 2 „ symmetry) and c-CC 13 CH 2 ( C s symmetry) have been 
observed in the 1-3 mm band toward the low-mass star-forming region L1527. 
We have detected 7, 3, and 6 lines of c-C 3 H 2 , c- 13 CCCH 2 , and c-CC 13 CH 2 , respec¬ 
tively, with the Nobeyama 45 m telescope, and 34, 6 , and 13 lines, respectively, 
with the IRAM 30 m telescope, where 7, 2, and 2 transitions, respectively, are 
observed with the both telescopes. With these data, we have evaluated the col¬ 
umn densities of the normal and 13 C isotopic species. The [c-C 3 H 2 ]/[c- 13 CCCH 2 ] 
ratio is determined to be 310 ± 80, while the [c-C 3 H 2 ]/[c-CC 13 CH 2 ] ratio is de¬ 
termined to be 61 ± 11. The [c-C 3 H 2 ]/[c- 13 CCCH 2 ] and [c-C 3 H 2 ]/[c-CC 13 CH 2 ] 
ratios expected from the elemental 12 C/ 13 C ratio are 60-70 and 30-35, respec¬ 
tively, where the latter takes into account the statistical factor of 2 for the two 
equivalent carbon atoms in c-C 3 H 2 . Hence, this observation further confirms the 
dilution of the 13 C species in carbon-chain molecules and their related molecules, 
which are thought to originate from the dilution of 13 C + in the gas-phase C + due 
to the isotope exchange reaction: 13 C + + CO —>■ 13 CO + C + . Moreover, the abun¬ 
dances of the two 13 C isotopic species are different from each other. The ratio 
of c- 13 CCCH 2 species relative to c-CC 13 CH 2 is determined to be 0.20 ± 0.05. If 
13 C were randomly substituted for the three carbon atoms, the [c- 13 CCCH 2 ]/[c- 
CC 13 CH 2 ] ratio would be 0.5. Hence, the observed ratio indicates that c-CC 13 CH 2 
exists more favorably. Possible origins of the different abundances are discussed. 
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Introduction 


Radio spectral lines of 13 C isotopic species of various molecules are now detectable 
with a reasonable observation time, thanks to recent developments of receiver and backend 
technologies. So far, spectral lines of 13 C isotopic species of abundant molecules, whose 
main isotopologue (normal species) lines are optically thick, have widely been observed to 
derive the column density of the normal species accurately. In this case, the 12 C/ 13 C ratio 
of each species is usually assumed t o be the same as the elemental 12 C/ 13 C ratio, which 


1998 


Milam et ah 


2005 ). However, 


is 60-70 in the solar neighborhood (Lucas & Liszt 
anomalies of the 12 C/ 13 C ratio for carbon-chain molecules have recently been recognized in 
nearby cold and dense molecular clouds. The anomalies include different 12 C/ 13 C ratios for 
different carbon atoms in a single molecular species and significantly high 12 C/ 13 C ratios 
for several carbon-chain molecules in comparison with the elemental 12 C/ 13 C ratio. 


As for the former anomaly, T akano et ah] (119981 1 reported the non-equivalent abundances 
of the three 13 C species of HC 3 N toward the cold dark cloud TMC-1 (Cyanopolyyne peak; 
CP). The abundance of the HCC 13 CN is higher than that of H 13 CCCN and HC 13 CCN 


by a factor of 1.4. 


Sakai etal. 


(120071 1 also reported a large difference between the 13 CCS 


and C 13 CS abundances toward TMC-1 (CP), the ratio of C 13 CS/ 13 CCS being 4.2. A 
similar non-equivalence of the_ 13 C isotop ologue s of a single molecular species is reported for 


CCH, C 3 S, and C 4 H (1 Sakai et al 


2010 a, 


2013). It is proposed that these non-equivalent 


abundances of the 13 C species reflect production pathways of each molecule and/or originate 

’nVionmnrr tVm 13/^ nrvoifmn unf Viin o m/vlo/^nlo 


from gas phase reac 

■ions 

2010 a; 

Furuva et al. 

2011 ) 


The second anomaly is that the 12 C/ 13 C ratios of some carbon-chain molecules are 
significantly higher than the elemental 12 C/ 13 C ratio (60-70) in the solar neighborhood, 
indicating that the 13 C isotope is diluted in these molecules. The most striking result is 
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reported for CCH (Sakai et ah 2010afh The CCH/ 13 CCH and CCH/C 13 CH ratios are higher 
than 250 and 170 toward TMC-1 (CP), respectively. High 12 C/ 13 C ratios are also seen in 
CCS, C^S, and CRT Such a dilution of the 13 C isotopic species was theoretically predicted 


by 


Langer et ah 


(1l984h . and it was confirmed for some carbon-chain molecules by the above 


observations. On the other hand, 


Liszt & Ziurvsl ( 20121 ) reported the 12 C/ 13 C ratios of CS, 


HNC, and H 2 CS to be in the range from 43 to 79 toward TMC-1 (CP), which is close to 
the elemental 12 C/ 13 C ratio. Hence, the degree of the 13 C dilution seems dependent on 
molecular species, although its origin is still controversial. 

It is therefore important to explore whether the 13 C isotope anomalies are seen in other 
molecules. Furthermore, the above observations were mostly carried out toward TMC-1 
(CP), although the anomaly was reported for CCS toward the cold dark cloud L1521E 
and for CCH toward the low-mass star-forming region L1527. With these in mind, we 
here report on the observations of the 13 C isotopic species of C-C 3 H 2 toward the low-mass 
star-forming region L1527. 

C-C 3 H 2 is a carbon-chain related molecule, a nd is known to exist widely in interstellar 


Madden et al. 

1989; 

Cox et al. 

1989) 


1989 b It has two kinds of 13 C isotopic species, 


clouds (e.g., 

as shown in Figure [0 The off-axis 13 C species (referred to hereafter as c-CC 13 CH 2 ) is 
statistically favored by a factor of 2 , because there are two equivalent carbon atoms in 


c-CsHofor 13 C subst 

itution for this species. It was detected in TMC-1 and other sources 

(Madden et al. 

1986; 

Gomez-Gonzalez et al. 

1986; 

Soezzano et al. 

2013 

). On the other 


hand, detection of the on-axis 13 C species (referred to hereafter as c- 13 CCCH 2 ) has not been 
r eported in cold_clouds as far as we know, although detection of one line was claimed by 


Gomez-Gonzalez et ab 


(119861) toward the Galactic center cloud Sgr B2. During the spectral 


line surveys of L1527 with Nobeyama 45 m telescope (a legacy program of Nobeyama Radio 
Observatory) and IRAM 30 m telescope (the IRAM large program ASAI), we detected 
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many spectral lines of the off-axis and on-axis 13 C isotopic species as well as those of the 
main isotopologue. On the basis of these extensive datasets, we have explored the 13 C 
abundance anomaly of C-C 3 H 2 in L1527. 


2. Observation 


2.1. Observation with Nobeyama 45 m Telescope 


The spectral line survey of the 3 mm band (80 GHz to 116 GHz) was carried out with 
Nobeyama 45 m telescope (referred to hereafter as NRO 45 m) toward L1527 from 2006 
to 2011, as the legacy project of the Nobeyama Radio Observatory. Transition frequencies 
and quantum numbers of the observed lines for C-C3H2 and its 13 C species are listed in 
Tables 1-3. The observed position was (« 2000 ) ^ 2000 ) = (04 h 39 m 53!89, 26°03 , ll , . / 0). In this 
observation, we used two SIS mixer receivers, S80 and S100, simultaneously in 2006-2008 


and the dual-polarization sideband-separating SIS receiver T100H/V f lNakaiima et ah 2008) 


in 2009-2011. The system temperatures varied from 250 to 350 K in 2006-2008 and from 
150 to 250 K in 2009-2011. The pointing of the telescope was checked every 1 ~ 1.5 hours by 
observing nearby SiO maser sources (NML-Tau, Ori-KL), and was maintained to be better 
than 6 ". The intensity scale was calibrated by the chopper-wheel method, and is estimated 
to be accurate within 20%. The main-beam efficiency is 0.43 in 2006-2008, 0.49 in 2009, 
0.42 in 2010-2011, as reported in the NRO website b The beam size is 19" at 86 GHz 1 . I 11 
2006-2011, we used a bank of acousto-optical radio spectrometers (AOSs) whose bandwidth, 
resolution, and channel spacing are 250 MHz, 250 kHz, and 125 kHz, respectively. In 2012, 
we used a bank of autocorrelators, SAM45, whose bandwidth, resolution, and channel 
spacing were set to be 500 MHz, 244 kHz, and 244 kHz, respectively. The frequency 


1 http://www.nro.nao.ac.jp/~nro45mrt/html/prop/status/Status_R14.html 
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resolution corresponds to ~ 0.8 km s^ 1 at 90 GHz. This resolution is not sufficient, because 
the typical line width toward this source is from 0.5 to 0.8 km s” 1 . However, we chose this 
resolution to cover the whole 3 mm band within the limited observation time. Although the 
lines are partly frequency-diluted, the integrated intensity is reliable. A position switching 
mode with the off-position, ( 02000 , £ 2000 ) = (04 h 42 m 35!9, 25°53'23 / /3), was employed in all 
the observations. 


2.2. Observation with IRAM 30 m Telescope 

Observations with IRAM 30 m telescope (referred to hereafter as IRAM 30 m) at 
Pico Veleta were performed in several runs between 2012 and 2013, as part of the ASAI 
(Astrochemical Surveys at IRAM) Large Program. The survey covers the spectral bands at 
3 mm (80-112 GHz), 2 mm (130-173 GHz), and 1 mm (200-276 GHz). The beam size is 
29", 17", and 12" at 86, 145, and 210 GHz, respectively, as reported in the IRAM websitt 2 . 
The wobbler switching mode was employed with a beam throw of 180". This amount of the 
beam throw is enough, because most of the emission comes from the 15"-20" diameter area 
around the pr otostar acco rding to the interferometer observations combined with single-dish 


observations (ISakai et ah 


2010blh The frontends were the broad-band EMIR receivers, and 


the backends were the FTS spectrometers in their 195 kHz resolution mode, corresponding 
to velocity resolutions of 0.21-0.73 km s -1 . The antenna temperature (T£) was converted 
to the main-beam brightness temperature (T m b) by dividing by B e g/F e g. Here, B e g and F e g 
are beam efficiency and forward efficiency, respectively. The B e g/F e g value is 0.85, 0.80, 
and 0.67 at 86, 145, and 210 GHz, respectively 2 . 


2 http: //www.iram.es/IRAMES / mainWiki/Iram30mEfficiencies 
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3. Results 

3.1. Overall results 

In the Nobeyama observations, we detected 7 lines of C-C 3 H 2 in the 3 mm band (Table 
[]J. Later, we further detected 8 , 9, and 17 lines in the 3 mm, 2 mm, and 1 mm bands in 
the ASAI observations. Examples of the observed spectra are shown in Figured! along with 
the best Gaussian fit to their profiles. In total, 41 lines of C-C3H2 were detected, among 
which 7 transitions were observed both in the Nobeyama and ASAI observations. Even 
high excitation lines of 8i i8 -8 0 , 8 and 8 2 ,7-8 i i8 , whose upper state energies are as high as 77 
K (53 cm -1 ), were detected. The detected lines are summarized in Tables |T] and [4] and are 
also shown in the energy level diagram (Figure [3]). Since C-C 3 H 2 has a pair of equivalent H 
nuclei, the rotational levels are classified into the ortho and para species. Radiative and 
collisional interconversion between the ortho and the para species is strongly forbidden, 
and hence, the ortho and the para species behave as different molecular species. In our 
observations, 19 and 18 lines are detected for the ortho and para species, respectively, while 
the remaining 4 lines detected in the 1 mm band are blended lines of the ortho and para 
species. 

In addition to the lines of normal species, we also detected the lines of the 13 C 
isotopic species. The line parameters of c-CC 13 CH 2 and c- 13 CCCH 2 species (Figured]) are 
summarized in Tables [2] and |3] respectively. For c-CC 13 CH 2 , we detected 6 lines in the 3 mm 
band in the Nobeyama observations, and 2, 4, and 7 lines in the 3 mm, 2 mm, 1 mm bands, 
respectively, in the ASAI observations. In total, we detected 19 lines for c-CC 13 CH 2 , as 
summarized in Tables [ 2 ] and [4] In contrast to the normal species, c-CC 13 CH 2 does not have 
the ortho and para classification, because the two hydrogen nuclei are no longer equivalent. 
For c- 13 CCCH 2 , we detected 3 lines in the 3 mm band in the Nobeyama observations, and 
2, 1, and 3 lines in the 3 mm, 2 mm, 1 mm bands, respectively, in the ASAI observations, 


as summarized in Tables [3] and [4] Hence, c- 13 CCCH 2 is definitively detected, as shown in 
Figure [21 This species has the ortho and para states as in the case of the normal species. 
It should be noted that c-CC 13 CH 2 is statistically favored by a factor of 2 in comparison 
with c- l 3 CCCH 2 , because two equivalent carbon atoms are available in C-C 3 H 2 for 13 C 
substitution for c-CC 13 CH 2 while only one for c- 13 CCCH 2 (Figured]). 


3.2. LTE analysis for c-C 3 H 2 


Assuming the local thermodynamic equilibrium (LTE) condition, we determine the 
rotation temperatures and the beam-averaged column densities of c-C 3 H 2 and its 13 C 
species by using the least-squares analysis with the following equations: 


m hv 
T h = — 


and 


r = 


k \_exp(hv/kT rot ) — 1 exp(hv / kT^ g ) — 1 

87 r 3 *V 


(1 


3hAvU(T mt ) |_ K kT, 


exp 


hv 


- rot 


- 1 


exp 


K 


kT r 


rot 


N, 


where T f, is the brightness temperature, h the Planck constant, v the transition frequency, 
k the Boltzmann constant, T rot the rotation temperature, 7b g the cosmic microwave 
background temperature, r the optical depth, S the line strength, p the dipole moment, 
Av the line width, U(T rot ) the partition function, E u the upper state energy, and N the 
total column density. The partition function is numerically calculated from energies and 


degener acies of the rotatio nal levels, which are taken from the spectr oscopy data 


CDMS ( 

Mull 

er et ah 

2005 

) (See 

Vrtilek et ah 

(1987 

); 

Boeev et ah 


so original spectr oscop y references: 


j ase, 


19871 ): 


Sp ezzano et ah 


( 2012 ).). 


Bogey et. ah 


(119861) : 


Because the size of the emitting region is comparable to or smaller than the telescope 
beam size, the effect of beam dilution has to be taken into account. This correction for the 
beam dilution is particularly important in this analysis, because the 1, 2 and 3 mm data, 
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which result in different beam sizes, are used simultaneously. Furthermore, the correction 
is also necessary for the 3 mm data taken with the two different telescopes. The correction 
for the beam dilution is made by dividing the observed intensity by the beam dilution 
factor 77 = dg/(0g + d'l), where 9 S and 9 b denote the source diameter and the antenna 
beam width (FWHM), respectively. This correction factor assumes a Gaussian beam and 
a Gaussian distribution of the source emission. The FWHM beam widths are assumed to 
be 1.7 x 10 3 /z/(GHz) and 2.4 x 10 3 /zz(GHz) in arcseconds for NRO 45 m and IRAM 30 nr 
respec tively, on the basis of the beam sizes mentioned in Sections 2.1 and 2 . 2 . 


Sakai et al. 


fl2010bi ) reported that C-C 3 H 2 is distributed within a diameter of 2000-3000 AU from the 


protostar, which corresponds to the size of 15 ,/ -20". Thus, the size of the emitting region 
is assumed to be 20". The intensities of the high excitation lines ( E a > 30 K) in the 3 
mm band observed with NRO 45 m are brighter by a factor of 1.3-2.4 than those observed 
with IRAM 30 m, which further supports such a compact distribution. In the least-squares 
analysis, the ortho and para species are treated as different molecules, and the ortho-to-para 
ratio is taken as a fitting parameter. This analysis is possible owing to the large dataset 
of the observed line intensities. Relative weights of the data are determined by the errors 
of the line parameters including the uncertainty of the intensity calibration (20%). The 
total (ortho plus para) column density is derived to be (1.08 ± 0.16) x 10 14 cm -2 . The 
ortho-to-para ratio is evaluated to be 3.1 ±0.4, which is comparabl e to t he stati stical weight 


of 3. This ratio is consistent with the previous report by 


Takakuwa et al. 


(1200lh (2.5 ±0.5). 


The rotation temperature is 


L1527 of 12.3 K ( Saka i et al. 


8.3 ± 0.3) K, which is lower than that reported for C 4 H 2 in 


2008). The errors denote the standard deviation of the fit. 


We also estimated the column densities of the two 13 C species, c- 13 CCCH 2 and 
c-CC 13 CH 2 , by assuming LTE condition. The column density and the rotation temperature 
of c- 13 CCCH 2 are evaluated to be (3.5 ± 0.7) x 10 11 cm -2 and ( 8.6 ± 1.5) K, respectively, 
by assuming the same ortho-to-para ratio as the normal species. On the other hand, 
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the column density and the excitation temperature of c-CC 13 CH 2 are determined to be 
(1.79±0.19) x 10 12 cm -2 and (7.3±0.4) K, respectively. The abundance ratio of [c-C 3 H 2 ]/[c- 
13 CCCH 2 ] is thus determined to be 310 ± 80, whereas that of [c-C 3 H 2 ]/[c-CC 13 CH 2 ] is 
61 ± 11. [c-CC 13 CH 2 ] represents the total abundance of the c-CC 13 CH 2 species without any 
correction for the statistical factor of 2 (Figured! Section 3.1). These results are higher 
than the ratios expected for the condition that 13 C is randomly distributed among the three 
carbon atoms and the elemental abundance ratio of 12 C/ 13 C is 60-70, as shown in Tabled) 

Figure d] shows the residuals in the least-squares fit divided by the uncertainty ST of 
each line, (T obs — T &t )/ST, where T obs and T fit are the observed intensity and the best-fit 
value, respectively. The uncertainty of the intensity ST for each line is estimated by taking 
into account the Gaussian fitting error of the line profile and the calibration uncertainty. 
Figure 0] shows a slight systematic residuals in the fit. Intensities of the lines with high and 
low upper state energies systematically show positive residuals, while those of the lines with 
moderate upper state energies tend to have negative residuals. Such a systematic trend 
can also be seen, even if the different source sizes of 15" and 30" are assumed. This may 
indicate that the source includes two components having different rotation temperatures. 
The physical structure of the source (i.e. the gradient in the density and the temperature) 
may also play a role. However, the systematic residuals are eliminated by a non-LTE 
analysis described below, and we therefore consider that the systematic residuals mainly 
originate from imperfection of the LTE assumption. 


3.3. LVG analysis for c-C 3 H 2 


We employ a lar ge ve loci ty g radient (LVG) analysis based on the non-LTE radiative 


transfer code Radex (van der Tak et al. 


20071) for th e analysis of the n ormal species. If 


the gas kinetic temperature is assumed to be 25 K flSakai et ah 


2010bl ). the H 2 density 
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and the total column density of C-C 3 H 2 are determined to be (2.9 ± 0.5) x 10 5 cm -3 and 
( 8.8 ± 1.1) x 10 13 cm -2 , respectively. Even if the gas kinetic temperature is assumed to 
be 10 K for reference, the total column density of C-C 3 H 2 does not change very much 
((1.18 ± 0.11) x 10 14 cm -2 ). Figure [5] shows the residuals of the fit for the gas kinetic 
temperature of 25 K, calculated in the same way as in Figure 01 The systematic errors seen 
in Figure Q] are now eliminated in Figure 0 Although the non-LTE analysis is favorable as 
seen here, the collisional cross sections are not reported for the 13 C isotopic species. Hence, 
we cannot determine the abundance ratios of the 13 C isotopic species relative to the normal 
species by using the non-LTE analysis. Since the total column density derived from the 
non-LTE analysis only differs from that derived from the LTE analysis by 20%, we employ 
the LTE results for discussions of the 12 C/ 13 C ratio of C-C3H2. 


4. Discussion 


4.1. Dilution of the 13 C species 


The 12 C/ 13 C ratios of C-C3H2 derived from the column densities are summarized in 
Table [5j The averaged 12 C/ 13 C ratio of C-C3H2 is calculated as 

3[c-C 3 H 2 ] 


Rav [c- 13 CCCH 2 ] + [c-CC 13 CH 2 ] 


( 1 ) 


The numerical factor 3 is introduced because three carbon atoms are available for the 13 C 
substitution. The R av values are evaluated to be 150 ± 30 under the LTE condition. As 
shown in Table 0 the 12 C/ 13 C ratios of C-C 3 H 2 are thus found to be higher than those 


expected from the interstellar elemental abundance ratio of 60-70 flLucas fe Liszt 


Milam et ah 


1998 


200 5). For comparison, we have also derived the 12 C/ 13 C ratio in L1527 by 


using the J = 1-0 line data of C 18 0 and 13 C 18 0 obtained in the ASAI and Nobeyama 
observations to be 68 ± 19. The R av value obtained above is also higher than this ratio. 
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In dense molecular clouds, the main reservoir of 13 C is 13 C0, which is the most 
abundant carbon-bearing molecule. As carbon-chain molecules are produced by reactions 
starting from C + , their 13 C species are formed by the reaction from 13 C + . Although the 
13 C + ion is produced from 13 CO with the reactions of He + , as 12 C + from 12 CO, 13 C + will go 


back to 13 CO by the following isotope exchange reaction ( Langer et ah 
2015 ): 


19841; 


R.oueff et ah 


13 1 


C + + 12 CO ->• ia CO + rz C + + 35 K 


13 


12 / 


( 2 ) 


This is the most important reaction for 13 C + , because the reaction with H 2 is a very slow 
radiative association. The electron recombination is also much slower than the reaction 
with 12 CO. This exchange reaction is exothermic, and hence, its backward reaction is slow 
in cold clouds. Then, the 12 C + / 13 C + ratio tends to be higher in cold clouds. The 13 C species 


of various molecules proc 
was initially proposed by 


uced from 13 C + thus become less abundant. This mechanism 


Langer et ah 


(119841 ). Woo ds fe Willac y! ( 20091 ) and 


Furuva et ah 


(20111) confirmed this trend by chemical model calculations. It should be noted that the 


isotope selective photodissociation, which could cause the 12 C/ 13 C isotope anomaly in 
photodissociation regions, is not important in our observations, because we are observing a 
dense part near the protostar. 

Such a dilution process can occur even in the protostellar envelope of L1527, whose gas 


kinetic temperature is 25 K (1 Sakai et ah 


2010 b, 


2 014 ). In order to verify this, we consider 


the following relation under the steady-state approximation: 

[ 13 C+] _ k h { 13 CO] 

[ 12 C+] “ k { [ 12 CO] ’ 

where k{ and kb are the rate coefficients of the forward and backward reactions of 


( 3 ) 


respectively, and [A"] denotes the abundance of A". These rate coefficients fulfill the following 
relation: 


kb = h exp — 


AG 

W 


( 4 ) 
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where AG represents the free energy difference between the right and left hand sides of 
the reaction (J2D- Since k^/kf is evaluated to be 0.2 at 25 K, the dilution process of 13 C + is 
barely possible in the warm region of L1527. 

L1527 is a warm-carbon-chain-chemistry source, where Cff 4 evaporated from dust 
grains in the warm (> 25 K) and the dense part near the protostar triggers efficient 
production of carbon-chain molecules including C-C 3 H 2 . CH4 is thought to be formed by 
hydrogenation from the neutral carbon atoms depleted onto dust grains in the cold starless 
core phase and/or the less dense phase. The 12 C/ 13 C ratio of the carbon atom is also 
expected to be larger than the elemental ratio as the 12 C + / 13 C + ratio, because the neutral 
carbon atom (C) is formed from electron recombination of the carbon ion (C + ). ffence, the 
12 C/ 13 C ratio of CH 4 formed on dust grains would also be higher. This may also contribute 
to the dilution of the 13 C species of c-C 3 H 2 . 


4.2. Abundance difference between two 13 C species 


As described in the Introduction, it is reported that several carbon-chain molecules 
show the abundance anomaly of the 13 C species. For CCH, CCS, C 3 S, C 4 H, and 
HC 3 N, the 13 C abundances for different carbon atoms in a single molecular species are 
different from one another. In order to compare the abundance of the two 13 C species 
(Section 3.1), we calculate the [c- 13 CCCH 2 ]/[c-CC 13 CH 2 ] ratio to be 0.20 ± 0.05. If 13 C 
were randomly substituted for the three carbon atoms, the [c- 13 CCCH 2 ]/[c-CC 13 CH 2 ] 
ratio would be 0.5. Hence, the observed ratio indicates that c-CC 13 CH 2 exists more 
favorably than c- 13 CCCH 2 . It has been thought that the main production pathway for 
c-C 3 H 2 is the electron recombination of C 3 Hjj\ This is because the cyclopropenyl cation 
C 3 H 3 + is a major ionic species, which is produced by the gas-phase reactions starting 


from ('11 1 (Sakai & Yamamoto 


20131 ). Since C 3 Hjj" has three equivalent carbon atoms, 
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C-C3H2 produced from C3H3 cannot contribute to the abundance anomaly of 13 C species. 
Hence, other reactions have to be considered as the dominant production pathway (e.g., 
C 2 H 2 + CH —> C 3 H 2 + H). 


Another possible reason for the abundance anomaly of the two 13 C species is the 
exchange reaction between c- 13 CCCH 2 and c-CC 13 CH 2 , although it is uncertain whether 
this reaction happens for the closed shell molecule, C-C 3 H 2 . The zero-point vibrational 
energies of c- 13 CCCH 2 and c-CC 13 CH 2 are evaluated to be lower than the normal_species by 
8.74 K and 53.2 K, respectively, by using the molecular constants by iDateo & Lee ( 19971 ). 


Hence, it would be possible that c-CC 13 CH 2 is enhanced relative to c- 13 CCCH 2 , if the 
exchange reaction is possible. In case of 13 CCH and C 13 CH, the reaction such as 


13 


CCH + H C 13 CH + H 


(5) 


is proposed flSakai et al. 


2010a). It is thus interesting to explore whether a similar process 


is possible for c-C 3 H 2 . It should be noted that the isotopologue with the lowest zero-point 


vibrational energy tends to 


lave higher abunc 


trend is uncertain for C 4 H (1 Sakai et al 


20071. 


ances 


2010 a. 


contribution of the position exchange reactions ( ISakai et al. 


'or CCH, CCS, and C 3 S, although this 
20131). This fa c t may suggest poss ible 


2010 a; 


Furuva et ah 


2011 ). 


More detailed studies of these reactions as well as chemical model calculations for the 13 C 
species in molecular clouds are awaited. 


5. Concluding Remarks 

In this study, we accurately determined the [c-C 3 H 2 ]/[c- 13 CCCH 2 ] and [c-C 3 H 2 ]/[c- 
CC 13 CH 2 ] ratios on the basis of the extensive data of the Nobeyama and ASAI observations. 
We found a significant dilution of the 13 C species, and a significant difference between the 
[c-C 3 H 2 ]/[c- 13 CCCH 2 ] and [c-C 3 H 2 ]/[c-CC 13 CH 2 ] ratios, even when the statistical factor 
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of 2 for the c-CC 13 CH 2 species is taken into account. This further confirms that the 13 C 
anomaly is an important phenomenon in astrochemistry. 


It should be noted that determination of the 12 C/ 13 C ratio requires accurate 
determination of the column density of the normal species. This was indeed possible in 
this study, where a large number of lines can be used to conduct a fine analysis. However, 
this is not always the case. On the other hand, the anomaly in the 12 C/ 13 C ratio among 
the different carbon atoms in a single molecule can be studied only by observations of the 
13 C species, which often emit optically thin lines. If the cause of the 12 C/ 13 C anomaly in 
a single molecule is established, it could be a new tracer to understand chemical processes 
and physical conditions of molecular clouds. 

The anomalies of the 13 C species found for c- CbH? and other carbon-chain molecules 
might be related to the 13 C anomaly in meteorites (F 


anomaly in dust forming regions of AGB stars (e.g., 


oss et al. 

2006 

). A 

'though the 12 C/ 13 C 

Milam et al. 

2009: 

Kahane et al. 

1992: 


Tomkin et al. 


19761 ) mainly contributes to the anomaly in meteorites, the 12 C/ 13 C anomaly 


in cold molecular clouds may also be responsible. Hence, understanding the abundance 
anomaly of the 13 C species and its evolution in star-forming regions is of particular 
importance to trace back from our planetary system to its interstellar origins. 
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Fig. 1.— Molecular geometry of C-C 3 H 2 and its 13 C species 
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Fig. 2.— Observed line profiles of C-C 3 H 2 and its 13 C species in L1527. The results of 
Gaussian fitting are also shown in red. 
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Fig. 3.— The energy level diagram of C-C 3 H 2 . 
by downward arrows. 


The transitions observed in L1527 are shown 
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Fig. 4.— Residuals of the LTE fit divided by uncertainty of each line, arranged in ascending 
order of the upper-state energy. 
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Transition 


Fig. 5.— The same as Figure [4j but obtained by assuming the LVG condition 





























Table 1. Observed line parameters of C-C 3 H 2 toward L1527 


Frequency 

(GHz) 

Transition 

o/p 

E u 

(K) 

<S M 2 ) 

(Debye 2 ) 

AV L SR a 
(km s -1 ) 

At> a 

(km s -1 ) 

rms 

(mK) 

T mb a 

(K) 

fT mb dv 
(K km s _1 ) 

CLTE b 

CLVG C 

telescope 

80.723180 

4 2 ,2-4i,3 

p 

28.8 

19.5 

-0.07(4) 

1.72(9) 

5.8 

0.12(3) 

0.25(5) 

0.48 

1.54 

NRO 






-0.09(6) 

0.81(10) 

5.9 

0.14(3) 

0.11(2) 

-0.30 

0.99 

IRAM 

82.093544 

2 0,2-ll,l 

p 

6.4 

14.6 

-0.15(13) 

1.40(3) 

34.3 

1.5(3) 

2.5(5) 

1.11 

0.43 

NRO 






-0.030(4) 

0.820(11) 

3.9 

1.8(4) 

1.6(3) 

2.09 

1.38 

IRAM 

82.966200 

3i,2-3o,3 

0 

16.0 

31.3 

-0.101(7) 

1.255(17) 

15.2 

1.1(2) 

1.6(3) 

0.56 

1.45 

NRO 






-0.063(4) 

0.788(5) 

3.5 

1.4(3) 

1.2(2) 

1.73 

2.16 

IRAM 

84.727688 

32,2-3l,3 

p 

16.1 

10.3 

0.150(5) 

1.154(12) 

5.5 

0.52(11) 

0.67(13) 

1.02 

1.17 

NRO 






-0.037(5) 

1.196(10) 

2.7 

0.42(8) 

0.54(11) 

1.90 

2.03 

IRAM 

85.338894 

2l,2-lo,l 

0 

6.4 

48.1 

0.340(6) 

1.239(14) 

21.1 

2.8(6) 

3.8(8) 

0.61 

-2.84 

NRO 






-0.070(5) 

1.065(16) 

2.9 

2.8(6) 

3.3(7) 

2.14 

-0.53 

IRAM 

85.656431 

4 3,2-42,3 

0 

29.1 

56.2 

-0.010(9) 

1.16(2) 

7.4 

0.46(9) 

0.57(11) 

-0.44 

2.94 

NRO 






-0.060(9) 

1.15(3) 

2.1 

0.26(5) 

0.33(7) 

-1.34 

2.60 

IRAM 

87.435318 

54,2-53,3 

p 

45.3 

26.2 

-0.2(2) 

2.0(5) 

10.5 

0.032(9) 

0.045(19) 

1.08 

1.73 

NRO 






0.00(11) 

1.9(3) 

1.8 

0.016(4) 

0.036(8) 

0.25 

1.42 

IRAM 

90.344082 

65,2-64,3 

0 

64.7 

98.3 

-0.16(19) 

1.2(5) 

3.8 

0.014(5) 

0.017(7) 

0.71 

0.81 

IRAM 

145.089606 

3l,2-22,1 

0 

16.0 

40.0 

0.004(2) 

0.727(5) 

12.8 

3.0(6) 

2.3(5) 

0.98 

0.30 

IRAM 

150.436555 

2 2,0- 1 l,l 

p 

9.7 

6.8 

-0.136(2) 

0.634(8) 

8.8 

1.7(3) 

1.2(2) 

0.080 

0.89 

IRAM 

150.820665 

4 0,4-3l,3 

p 

19.3 

36.9 

-0.162(4) 

0.779(10) 

8.2 

2.3(5) 

1.9(4) 

0.68 

-0.52 

IRAM 

150.851908 

4 i ,4 —3o,3 

0 

19.3 

110.7 

-0.201(13) 

0.92(3) 

8.3 

3.3(7) 

3.2(6) 

1.02 

-0.92 

IRAM 

150.954691 

62,4-61,5 

p 

54.7 

20.3 

-0.23(8) 

0.6(3) 

9.3 

0.05(2) 

0.025(9) 

0.31 

-0.63 

IRAM 

151.039173 

63,4-62,5 

0 

54.7 

60.8 

-0.10(4) 

0.71(9) 

8.8 

0.10(2) 

0.085(20) 

-0.010 

-0.70 

IRAM 

151.343878 

5l,4 _ 5o,5 

0 

35.4 

32.4 

-0.065(8) 

0.811(17) 

6.7 

0.43(9) 

0.37(8) 

-0.15 

0.79 

IRAM 

151.361105 

52,4“5i ,5 

p 

35.4 

10.8 

0.00(2) 

0.75(5) 

6.4 

0.14(3) 

0.11(2) 

-0.94 

-0.65 

IRAM 

155.518303 

32,2-2l,l 

p 

16.1 

17.8 

0.005(2) 

0.687(5) 

8.2 

2.0(4) 

1.5(3) 

0.38 

0.29 

IRAM 

204.788926 

42,2-33,1 

p 

28.8 

11.2 

-0.036(7) 

0.761(16) 

7.4 

0.52(10) 

0.44(9) 

-0.66 

-0.054 

IRAM 

216.278756 

33,0— 2 2,1 

0 

19.5 

45.6 

-0.128(7) 

0.892(15) 

5.2 

2.0(4) 

1.9(4) 

-1.71 

0.74 

IRAM 

217.822148 

60,6—51,5> 61,6-50,5 

p/o 

38.6 

58.3 

-0.092(3) 

0.922(8) 

4.4 

2.5(5) 

2.5(5) 

-0.42 

-2.22 

IRAM 

217.940046 

5i ,4—42,3 

0 

35.4 

110.3 

-0.057(2) 

0.863(5) 

5.5 

1.8(4) 

1.7(3) 

-0.25 

-0.073 

IRAM 

218.160456 

52,4~4l,3 

p 

35.4 

36.8 

-0.022(4) 

0.815(11) 

6.5 

0.73(15) 

0.63(13) 

-0.81 

-1.01 

IRAM 

218.732732 

72,6-7i,7i 7i,6-7o,7 

p/o 

61.2 

11.0 

0.05(3) 

0.92(8) 

4.4 

0.071(15) 

0.084(18) 

1.31 

-1.46 

IRAM 

227.169138 

4 3,2-32,1 

0 

29.1 

61.6 

-0.093(4) 

0.827(9) 

8.4 

1.8(4) 

1.6(3) 

-0.95 

0.96 

IRAM 

244.222150 

32,1 _ 2i,2 

0 

18.2 

7.3 

-0.047(3) 

0.761(6) 

6.0 

1.5(3) 

1.2(2) 

0.24 

1.40 

IRAM 

249.054368 

52,3-43,2 

0 

41.0 

76.3 

-0.050(4) 

0.791(10) 

5.1 

1.0(2) 

0.89(18) 

-0.18 

0.85 

IRAM 

251.314367 

71,7—60,6 > 7q,7-6i,6 

p/o 

50.7 

69.0 

-0.014(3) 

0.830(8) 

5.7 

1.6(3) 

1.5(3) 

1.19 

-1.86 

IRAM 



Table 1—Continued 


Frequency 

(GHz) 

Transition 

o/p 

E u 

(K) 

<s m 2 > 

(Debye 2 ) 

AVLSR a 
(km s -1 ) 

Ai> a 

(km s -1 ) 

rms 

(mK) 

T mb a 

(K) 

fTmbdv 
(K km s -1 ) 

<7LTE b 

CT LVG C 

telescope 

251.508709 

6 1,5 _5 2,4 

p 

47.5 

47.4 

0.029(9) 

0.74(2) 

6.0 

0.41(8) 

0.35(7) 

0.41 

-1.20 

IRAM 

251.527311 

6 2,5 _5 1,4 

0 

47.5 

142.2 

-0.005(4) 

0.791(10) 

6.4 

1.1(2) 

1.0(2) 

0.81 

0.15 

IRAM 

252.409837 

8 l,7 -8 0,8> 8 2,7 —8 1,8 

p/o 

76.5 

11.1 

-0.05(6) 

0.94(14) 

5.3 

0.051(12) 

0.038(10) 

3.41 

1.22 

IRAM 

254.987652 

5 3,3 _4 2,2 

p 

41.1 

26.8 

0.048(8) 

0.752(18) 

5.6 

0.36(7) 

0.30(6) 

-1.25 

-1.04 

IRAM 

260.479746 

53,2-44,1 

0 

44.7 

25.8 

0.00(11) 

0.83(3) 

6.8 

0.32(6) 

0.31(6) 

0.62 

0.68 

IRAM 

261.831811 

33,l- 2 0,2 

p 

19.0 

1.4 

0.02(3) 

0.68(6) 

18.4 

0.30(6) 

0.19(4) 

-1.38 

0.91 

IRAM 

265.759481 

44,1-33,0 

0 

32.2 

89.8 

-0.076(5) 

0.879(13) 

8.2 

1.4(3) 

1.4(3) 

-2.33 

0.80 

IRAM 


Note. — The numbers in parentheses represent the errors in unit of last significant digits. 
a Obtained by the Gaussian fit. 


’The residual divided by the intensity error, shown in Figure f4l 
'The same as <jltE> but shown in Figure [ 5 ] 


to 

Or 






Table 2. Observed line parameters of c-CC 13 CH 2 toward L1527 


Frequency 

(GHz) 

Transition 

E u 

(K) 

<s M 2 ) 

(Debye 2 ) 

AVLSR a 
(km s -1 ) 

Av a 

(km s _1 ) 

rms 

(mK) 

7 mb 

(K) 

fT^dv 
(K km s -1 ) 

_ b 

CT LTE 

telescope 

80.047537 

2o,2~1i,1 

6.3 

27.8 

-0.07(2) 

1.25(5) 

7.5 

0.16(3) 

0.23(5) 

1.68 

NRO 

80.775347 

3i,2-3o,3 

15.7 

21.0 

0.19(16) 

1.5(4) 

6.4 

0.021(6) 

0.036(11) 

0.29 

NRO 

83.474137 

32,2-3l,3 

15.9 

20.5 

-0.00(16) 

1.5(4) 

6.8 

0.022(7) 

0.033(12) 

0.35 

NRO 





0.00(2) 

1.1(5) 

3.5 

0.018(5) 

0.027(7) 

0.053 

IRAM 

84.185635 

21,2-3-0,1 

6.3 

31.3 

0.015(10) 

1.18(2) 

4.3 

0.21(4) 

0.28(6) 

1.65 

NRO 





-0.08(2) 

0.78(5) 

3.6 

0.26(5) 

0.22(4) 

2.29 

IRAM 

114.897371 

3o,3-2l,2 

11.8 

50.5 

0.00(4) 

0.90(6) 

19.6 

0.32(7) 

0.30(6) 

0.80 

NRO 

115.524356 

3l,3-2o,2 

11.8 

50.7 

-0.00(6) 

0.78(7) 

22.3 

0.31(7) 

0.32(7) 

-0.052 

NRO 

140.432881 

3l,2-22,1 

15.7 

23.9 

-0.13(3) 

0.53(10) 

7.3 

0.15(4) 

0.081(18) 

0.25 

IRAM 

147.702239 

2 2,0- 1 l,l 

9.6 

13.9 

-0.18(3) 

0.65(5) 

9.4 

0.17(4) 

0.10(2) 

0.021 

IRAM 

148.114191 

4l ,4 _ 3(3,3 

19.0 

71.6 

-0.161(18) 

0.78(4) 

6.9 

0.18(4) 

0.14(3) 

-0.40 

IRAM 

153.894626 

32,2-2l,l 

15.9 

34.8 

-0.10(2) 

0.63(6) 

9.2 

0.16(4) 

0.13(3) 

-0.83 

IRAM 

212.457672 

33,0-22,1 

19.2 

31.7 

-0.02(5) 

0.94(11) 

7.3 

0.081(18) 

0.083(19) 

-3.06 

IRAM 

213.843018 

6 1,6-5o, 5> 6 0,6- 5 l,5 

37.9 

113.5 

0.02(3) 

0.96(6) 

4.4 

0.088(18) 

0.085(18) 

0.66 

IRAM 

213.872779 

5l,4- 4 2,3 

34.8 

70.7 

0.06(5) 

0.79(11) 

4.4 

0.045(10) 

0.029(7) 

0.11 

IRAM 

214.313034 

52,4-4l,3 

34.8 

70.8 

-0.07(5) 

0.82(11) 

4.2 

0.052(12) 

0.060(13) 

0.78 

IRAM 

225.435934 

4 3,2-32,1 

28.6 

40.3 

0.02(9) 

1.00(2) 

7.9 

0.051(14) 

0.044(14) 

-0.027 

IRAM 

237.998109 

32,l-2l,2 

17.8 

5.1 

-0.08(6) 

0.71(14) 

4.3 

0.037(10) 

0.019(6) 

-0.53 

IRAM 

246.723236 

7 l,7- 6 0,6) 7 0,7-6l,6 

49.7 

134.3 

0.21(6) 

0.70(14) 

5.6 

0.044(12) 

0.040(10) 

1.25 

IRAM 


Note. — The numbers in parentheses represent the errors in unit of last significant digits. 
a Obtained by the Gaussian fit. 

k Residuals of the LTE fit divided by the uncertainty of each line. 



Table 3. Observed line parameters of c- 13 CCCH 2 toward L1527 


Frequency 

(GHz) 

Transition 

o/p 

E u 

(K) 

<Sp 2 > 

(Debye 2 ) 

AVLSR a 
(km s -1 ) 

Au a 

(km s -1 ) 

rms 

(mK) 

7 mb 

(K) 

fTmbdv 
(K km s -1 ) 

CLTE b 

telescope 

81.150881 

2 0 , 2 -ll,l 

p 

6.3 

31.0 

0.10(18) 

1.3(4) 

5.5 

0.016(5) 

0 . 021 ( 8 ) 

0.84 

NRO 






-0.00(3) 

0.9(3) 

4.1 

0.018(6) 

0.026(8) 

1.17 

IRAM 

82.303747 

2 1 , 2 - 1 o ,1 

0 

6.3 

96.2 

-0.17(12) 

1.1(3) 

12.6 

0.048(14) 

0.049(16) 

0.22 

NRO 






-0.08(7) 

0.83(8) 

3.9 

0.064(14) 

0.060(13) 

2.07 

IRAM 

114.381212 

3<),3 - 2 l,2 

0 

11.8 

159.0 

-0.15(11) 

1 . 1 ( 2 ) 

14.8 

0.08(2) 

0.07(3) 

0.30 

NRO 

145.353129 

3l,2“ 2 2,l 

0 

15.7 

96.1 

0 . 10 ( 11 ) 

0 . 6 ( 2 ) 

10.3 

0.044(16) 

0.033(13) 

-0.80 

IRAM 

211.613800 

33,0- 2 2,l 

0 

18.9 

75.0 

0.02(7) 

0.69(17) 

4.9 

0.033(10) 

0.024(7) 

-2.15 

IRAM 

212.387110 

6 0,6 _5 1,5> 6 1,6 _5 0,5 

p/o 

37.6 

117.2 

0.18(9) 

0 . 8 ( 2 ) 

6.5 

0.039(12) 

0.026(8) 

0.39 

IRAM 

245.044421 

7 1 , 7 - 60,65 7 0,7- 6 l,6 

p/o 

49.4 

138.6 

-0.04(10) 

1 . 2 ( 2 ) 

4.8 

0.029(8) 

0.033(10) 

2.72 

IRAM 


Note. — The numbers in parentheses represent the errors in unit of last significant digits. 
a Obtained by the Gaussian fit. 

k Residuals of the LTE fit divided by the uncertainty of each line. 



Table 4: Number of detected lines. 



NRO 45 m 

IRAM 30 

m 


Species 

3 mm 

3 mm 

2 mm 

1 mm 

Total 

c-C 3 H 2 

7 a 

8 

9 

17 

41 

c-CC 13 CH 2 

6 

2 b 

4 

7 

19 

c - 13 CCCH 2 

3 

2 b 

1 

3 

9 


“All the transitions are also detected with IRAM 30 nr. 
b All the transitions are also detected with NRO 45 nr. 
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Table 5: The observed 12 C/ 13 C ratios of c-C 3 H 2 . 


Ratio 

Observed 3, 

Expected 13 

[c-C 3 H 2 ]/[c-CC 13 CH 2 ] 

61 ± 11 

30-35 

[c-C 3 H 2 ]/[c- 13 CCCH 2 ] 

310 ± 80 

60-70 

3[c-C 3 H 2 ]/([c-CC 13 CH 2 ]+[c- 13 CCCH 2 ]) 

150 ± 30 

60-70 

[c- 13 CCCH 2 ]/[c-CC 13 CH 2 ] 

0.20 ±0.05 

0.5 


Note. - The errors represent the standard deviation. 

“Ratios directly derived from the column densities of C-C 3 H 2 , c- 13 CCCH 2 , and c-CC 13 CH 2 . 

^Ratios expected for the condition that 13 C is randomly distributed among the three carbon atoms and the 
elemental abundance ratio of 12 C/ 13 C is 60-70. 



